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Abstract 

We provide universal formulae for the limiting stretching and bend- 
ing energies of triangulated membrane networks endowed with near- 
est neighbor bond potentials and cosine-type dihedral angle poten- 
tials. The given formulae account for finite elasticity and solve some 
deficiencies of earlier results for Helfrich-type bending energies, due 
to shape-dependence and sensitivity to mesh distortion effects of the 
limiting elastic coefficients. We also provide the entire set of the elas- 
tic coefficients characterizing the limiting response of the examined 
networks, accounting for full bending-stretching coupling. We illus- 
trate the effectiveness of the proposed formulae by way of example, 
on examining the special cases of cylindrical and spherical networks 
covered with equilateral triangles, and discussing possible strategies 
for the experimental characterization of selected elastic moduli. 

Key words: Membrane networks, Dihedral angle potentials, Nearest 
neighbor bond potentials, Multiscale modeling, Limiting energies, Bending 
stiffness. 



1 Introduction 

Polyhedral membrane models endowed with bond and dihedral angle poten- 
tials (hereafter referred to as membrane networks) are widely employed in 
many different physical and engineering problems, including molecular dy- 
namics (MD) simulations of bio- and nano-structures (refer, e.g., to the re- 
cent review papers Midler et al., 2006, Tu et al, 2008 and references therein), 
and non- conforming models of plates and shells (Davini and Pitacco, 2000, 
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Angelillo et al. 2006), among others. Referring to the special case of the 
red blood cell membrane, the bond potentials reproduce the protein fila- 
ments forming the membrane skeleton, while the dihedral angle potentials 
mimic the bending rigidity endowed by the lipid bilayer (Tu ct al., 2008; 
Discher et al., 1998; Marcelli ct al., 2005; Dao and Suresh, 2006; Hale et al., 
2009; Hartmann, 2010). In such a case, the discrete nature of the net- 
work is motivated by the finite size of the protein filaments, and the con- 
venience of numerical approaches like MD and Monte Carlo simulations 
(Kroll and Gomppcr, 1992). On the other hand, the use of membrane net- 
works for the development of lumped strain models of plates and shells allows 
for capturing local effects due, e.g., to strain localization, plastic hinges and 
folding-hke buckling (Davini and Pitacco, 2000; Angelillo et al., 2006). 

A key question for such models is the determination of the continuum 
limit of stretching and bending energies, when the nearest neighbor distance r 
becomes infinitesimally small (Seung and Nelson, 1988; Gompper and Kroll, 
1996; Discher et al., 1997; Zhou and Joos, 1997; Schmidt, 2006; Schmidt, 
2008; Hartmann, 2010). While the limit form of the membrane energy is well 
established in the literature, at least for harmonic bonds (Scung and Nelson, 
1988; Discher et al., 1997; Zhou and ,Jo6s, 1997) some questions still remain 
open regarding the continuum limit of the bending energy, when cosine- 
type dihedral angle potentials are employed. As a matter of fact, ear- 
lier results in terms of Helfrich-type limiting energies for triangulated net- 
works suffer from deficiencies due to shape-dependence and mesh distortion 
effects (Scung and Nelson, 1988; Gomppcr and Kroll, 1996; Nelson ct al., 
2004; Espriu, 1987; Baillic at al., 1990). Different is the case of bending 
energies based on nodal averages of dihedral angles, which exhibit a more 
uniform limit behavior (Gompper and Kroll, 1996). 

We deal in the present paper with the derivation of limiting stretching 
and bending energies of static membrane networks undergoing large in-plane 
and out-of-plane deformations. We consider a discrete mass spring system 
with nearest neighbor and dihedral spring interaction subject to deforma- 
tion, and provide explicit formulae for the energy in the system when the 
number of particles (vertices) tends to infinity (i.e., r — 0). The limit- 
ing energy only depends on the geometric data of the deformation map- 
ping. As compared to previous works on the limiting bending energies of 
membrane networks (Seung and Nelson, 1988; Gompper and Kroll, 1996), 
the present study follows a different from discrete to continuum approach, 
by interpreting the discrete deformation of the network as the restriction 
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of a continuous deformation mapping / to a triangular lattice. This is in 
agreement with the Cauchy-Born rule of crystalline matter (Ericksen, 2008), 
and avoids mesh- dependent effects (Miillcr ct al., 2006; Gomppcr and Kroll, 
1996; Nelson ct al., 2004). In the bending-dominant deformation regime, 
when the deformed configuration of the network approximatively coincides 
with an equilateral triangulation of a curved surface, the formula we obtain 
for the limiting bending energy reduces to 
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where: D is the dihedral constant; U is the reference configuration: H is 
twice the mean curvature, and K is the Gaussian curvature. The formula 
above corresponds to a Helfrich-type energy J^{k,hH'^ + kgK) du, with 
Kg = —Akh/^ and kh = \/3-D/2 for any geometry, kh and kq being the 
bending and Gaussian rigidities, respectively (Hclfrich, 1973). This formula 
in particular applies to lipid membranes which show a fluid-like behavior 
in the in plane coordinates. Such a behavior leads to the relaxation of all 
stretching energy contributions and thus results in a pure bending energy 
functional. Those membranes can be described accurately within our frame- 
work of spring networks if we allow for a distribution of defects in the under- 
lying lattice network on a mesoscopic scale. In this way one can find meshes 
with almost no stretching contributions in any deformed configuration. Also 
note that the energy then only depends on the curvature tensor of the surface 
f{U) and / is merely some parameterization of that surface used to express 
the bending energy If{u) ~ dS in local coordinates. 

For general deformations of triangulated networks with both stretching 
and bending contributions, however, the limiting energy expressions will nat- 
urally depend on the lattice reference configuration. As a matter of fact, the 
prediction of mesoscopic values of the referential moduli from MD simula- 
tions, eventually different from region to region of the model, allows one to ad- 
equately inform simulations at the continuum scale, such as, e.g., Lagrangian 
finite element models. Moreover, stability analyses are naturally performed 
through referential (or thermodynamic) elastic constants, as shown, e.g., in 
Zhou and Joos, 1996; Zhou and Joos, 1997; Dischcr ct al., 1997; Hess ct al., 
1997. 

Concerning the present estimates of the bending and Gaussian rigidities, 
we wish to remark that Seung and Nelson (1988) found the different result 
Kg = —i^H and obtained a limiting bending energy of the form kh/2 j^{H'^ — 
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2K) du. The same authors also estimated = \/3D/2 by applying this 
formula to a cylindrical network. Subsequently, Gompper and Kroll (1996) 
made use of the result kg = —i^h to derive a different value of kh for the 
sphere. These authors computed the limiting value of a discrete notion of 
the Laplacian part of the Helfrich energy {(3^}"^^ = Jjj H'^ /2 du), and sub- 
tracted the total curvature of the sphere = K du = 4tt to fH-O"'^^. On 
matching kh^PI-O"''^^ — i^*"*) to the limiting value of the total dihedral energy 
of the sphere (4D7r/v^, cf. formula (19) of Gompper and Kroll, 1996), they 
finally obtained kh = a/S-D/S for such a geometry. 

It is easily seen that the present limiting bending energy J^^SH"^ — 
8K) du exactly reduces to the limiting dihedral energies of both the sphere 
and the cylinder, with kh = \/3D/2 in each case (cf. Sec. 6). We wish to 
remark that the scale bridging approach behind the 'universal' results = 
— 4k/^/3 and kjj = \piDj2 assumes that the discretized surface matches the 
limiting shape in correspondence with the vertices of the triangulation. In 
more detail, such results are obtained looking at the discrete dihedral en- 
ergy under the above Cauchy-Born type assumption, without making any 'a 
priori' guess on the kqI^^h ratio. Interestingly, the present estimates of kh 
and Kg correspond to those obtained in Lidmar et al., 2003 for the bending- 
dominant regime, through direct arguments. These authors 'assume' that 
the limiting bending energy can be written as a Helffrich-type energy with 
two free parameters (k// and kg, namely) in such a regime. Then they use 
the two explicit cases of a sphere and a cylinder to determine the above 
parameters. An analogous result is deduced in Fcdosov ct al., 2009 for the 
special case of a spherical shell. As we already observed, the present result 
for the bending energy is obtained without any a-priori assumptions on the 
structure of the limiting energy and the shape of the network and, as far as 
we know, represents the first complete proof to-date that the limiting contin- 
uum energy of membrane networks endowed with cosine-type dihedral angle 
potentials is actually a Helffrich type functional in the bending-dominant 
regime. 

The general formula we provide for the bending energy accounts for the 
coupling of large stretching and bending deformations, which is a key feature 
of finite deformation shell theories (Naglidi, 1972). It allows us to handle ar- 
bitrarily distorted triangulations, thus resolving the distortion sensitivity ef- 
fects that derive from the use of a pure bending energy in the continuum limit 
(Gompper and Kroll, 1996; Nelson et al., 2004; Espriu, 1987; Baillie at al.. 
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1990). Concerning the limiting stretching energy, our findings confirm the 
results given in Seung and Nelson, 1988; Discher et al., 1997; Zhou and Joos, 
1997. 

The paper is organized as follows. We first introduce the examined net- 
work model and the basic assumptions of the adopted scale bridging approach 
in Sec. 2. Next, we develop universal formulae for the limiting stretching and 
bending energies of such a network in Sees. 3 and 4, respectively. Elastic mod- 
uli are determined for the resulting effective theory in Sec. 5. In Sec. 6, we 
examine the special cases of cylindrical and spherical networks covered with 
equilateral triangles and undergoing pure bending. Sec. 7 discusses a mi- 
cropipette aspiration test as a nontrivial example for a general stretched and 
bent surface. We end by summarizing the main conclusions of the present 
work in Sec. 8. 



2 The model 



Let £ denote the triangular lattice C = [ ^ /^^^ ) Z in the Euclidean plane. 



1 1/2 \ 



V3/2) 



We consider a rescaled piece Cr = rCnU of C, where r > and f/ C open, 
which we consider as our reference configuration. A mapping f : U — )■ 
(which we assume to be sufficiently smooth up to the boundary) defines a 
discrete deformation when restricted to Cr (Fig. 1): 

Cr 3 U ^ f{u). 

The energy of the system shall be due to stretching and bending contri- 
butions: 

For large deformations we assume that the stretching energy is given by 
nearest neighbor interactions 

Ef-'^\f)= J2 Wr{\f{u)-f{v)\) (1) 

\u — v\ —r 

with Wr{s) = r'^W{s/r) for some smooth W. We note that such a scaling 
assumption is necessary in order to account for uniform response to stretching 
on the nearest neighbor scale |n — ?;| = r that is comparable to the size of 
an elementary lattice cell. For quadratic W{s) = ^{s — 1)^ one has Wr{s) = 
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Figure 1: Reference and deformed configurations of a membrane network 
generated by a triangular lattice. 



|(s — r)^. By way of contrast, in the regime of infinitesimal elasticity, where 
{Df)^Df — Id = 0(r), we rescale by setting Wr{s) = W{s/r) in order to 
obtain a finite energy density per unit area. 

The bending terms shall be given in terms of linear "cosine-type dihedral 
springs" , which associate energy to the cosine of the dihedral angle between 
two adjacent lattice cells. More precisely, we set 



^bend^^) = J2 (l - COS AA') 



A,A'GCr 
neighbors 



- E 

2 ^ 



(2) 



A,A'GCr 
neighbors 



where Cr denotes the set of equilateral triangles A C f/ of sidelength r with 
vertices in Cr', AA' is the dihedral angle between /(A) and /(A'); riA is the 
unit normal to /(A); and the summations run over all the pairs A, A' G Cr 
which share a common side. (Note that we tacitly assume that all /(A) are 
non-degenerate, the derivative Df has rank 2 everywhere and the normals 
are chosen with consistent orientation.) 

We will now consider the limiting energy expressions as r — separately. 
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3 Limiting stretching energy 



The computation of the continuum hmit r — of an energy of the form (1) 
is not new and rather straightforward. We include the short argument for 
the sake of completeness. 



3.1 Finite elasticity regime 

Denote by (-D/)a the average value of Df on A G and introduce the 
lattice vectors 

l/2\ / -1/2 \ /I 



Then (1) can be rewritten as 

AeCr 1 = 1 

where the error term 0{r) accounts for both the linearization error and possi- 
ble neglection of springs near the boundary. Since | A| = -\/3r^/4, a Riemann 
sum argument now gives 



^stretch. _ limEf^t^^r/) 

r— 5>0 

3 



-^Y. f W{\Df{u)a,\)du 



or, in terms of the first fundamental form = [D f (u))^ D f (u) , i.e., the 
right Cauchy-Green strain tensor, of the surface 

V3 .^^ Ju 

(■, ■) denoting the Euclidean scalar product. In coordinates = (gij) this 
expression becomes 



.stretch/ /.N 2 f f..rf V gu - 2^35-12 + Sg22 



= — [W 
v3 Ju 



W[ An + 2^3^12 + 3^22 , ^ ^^^^ , 
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Note that our derivation gives the exphcit error bound 

I ^stretch (y.) - ^«t>^^teh^y^| _ Q^^^ 

and it is not hard to see that this bound is optimal. 

3.2 Infinitesimal elasticity regime 

Here the rescaled Green- St. Venant strain tensor = '''^'^^ 2r^~^'^ 
der one and so / is an almost isometric immersion. More precisely, we as- 
sume that f = fr are lattice defomations such that E = lim^-i^oEj. exists. 
The contributions W{\{Df)/:^ai\) to the energy can be Taylor expanded as 
l^(|p/)Aa,|) = l^(l + r(a„(E,)Aa,)+0(r2)) = il^"(l)r2(a„ (E,)Aa,)2 + 
O(r^) and we obtain from (1) that 

3 



l^Ef^tc^ifr) = hmi Yl (E.)Aa.)^ + 0(r) 

A&Cr i=l 

= [ /iTr(E2) + -(TrE)2dM 
Ju 2 



is isotropic with Lame constants fi = \ = Such a result is in agree- 

ment with the predictions given in Seung and Nelson, 1988; Discher et al., 
1997; Zhou and Joos, 1997. Here the optimal error bound follows as 0{r + 
|E,-E|). 



4 Limiting bending energy 

Our derivation of the limiting bending energy is at variance with earlier com- 
putations (see Seung and Nelson, 1988; Gompper and Kroll, 1996; Discher et al., 
1997), and relies on the first and second fundamental forms of f{U). We de- 
fine the nearest neighbor vectors of the dual lattice as 

^^=(1/273)' ^^=(1/2^)' ^'=[-l/J)- 
Rewriting (2) as 

neighbor of A 
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we are led to compute the individual contributions \n/x — t^a'P for a fixed 
triangle A. After a rigid transformation, which in particular leaves the first 
fundamental form = [D f {u))^ D f {u) and the second fundamental form 



i=l 



\n[u] 



\duj{u)xdu2f(u)\ 



^ the unit normal) invariant, we may assume that the 



barycenter of A is and / = (/i, /2, /s) is such that /(O) = and -D/3 = 0. 
(See, e.g., Kuchncl, 2002 for a basic introduction to the differential geometry 
of surfaces.) 

After a tedious but straightforward computation, using that 

n(0) = (0,0,1) and so 
11(0) = ^^3(0), 

one finds 



\n/\ - riA' 



A'eCr 
neighbor of A 



2 

= ^^Lm (12(b., 11(0) b,) - Tr 11(0))^ + 0{r') 

A Riemann sum argument finally yields 



r— !>0 



U 



D A (a„Ia,)(12(b„IIb,)-TrII)' ^^ 

12^3^ detl 

1=1 



which in coordinates I(m) = {gij) and II(m) = {hij) can also be written as 



4V3det(^ij) (3) 
8gi2hiihi2 + 25(22(/iii + 3/1^2)) 
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Again our derivation yields the optimal error estimates 



The bending energy density is quadratic in II and can be expressed as 

in terms of the fourth order tensor © = (dijki) with dijki = dkiij, dijki = djM 
and 

'^1111 = 9ll + 2(^22; C^1212 = 9ll + 3fi'22) C?2222 = '^9ll^ 

dui2 = — 2(7i2, (^1122 = ~9U^ '^1222 = 0. 

This in turn reduces to (h, Dh) with h = (/in, /i22! 2/112) and the second 
order 3x3 tensor 

J-) /9ii + 25-22 -511 

2V3det{gij) y _2g^^ ° ' 

from which the nominal bending stresses and the the referential bending 
moduli can be readily read off by differentiation with respect to h. Note 
carefully that these will depend on the first fundamental form (5.^) thus 
introducing a coupling between stretching and bending contributions. Sym- 
metrically, the differentiation of the bending energy density with respect to 
Qij provides contributions to the nominal stretching stresses and referential 
stretching moduli. We remark that (3) accounts for possible heavy stretch- 
ing and shearing deformations of the reference lattice, due to the presence 
of the first fundamental form l{u) in the limiting energy. We devote the 
following section to the derivation of the complete set of elastic moduli that 
characterize the present network model. 

For isometric immersions, when I(m) = Id, formula (3) reduces to 

= ^ l^iSH' - 8K) du, (5) 

where H = TtII is twice the mean curvature and K = det II is the Gauss 
curvature of the surface, which provides a good approximation to the energy 
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also when I Id {bending- dominated deformation regime). It has to be 
noted that I = Id would imply that / is an isometric immersion and thus 
K = 0. However, general surfaces can be triangulated with I ^ Id if one 
allows for defects in the reference triangular lattice on a mesoscopic scale. 
More precisely: Choosing a partition of the surface into regions of area and 
boundary length scaling with resp. f for r ^ f ^ 1 one can parameterize 
the surface with almost trivial first fundamental form if one allows for defects 
in the reference lattice along the mesoscopic boundary of these elements. The 
number of atomic bonds across this boundary is 0((fr)^^), which is much less 
than the total number of bonds, which scales with r^^. In the limit r, f — )■ 
the contribution of the defects becomes negligible. Also the stretching terms 
will vanish as I approaches Id, and so the limiting continuum energy is given 
by (5). More precisely, the extra stretching energy from the lattice defects 
is bounded by 0(0, while |I — Id| = 0(r), leading to O(f^) contributions 
to ^''^ti'stch g^j^i^^ general, 0(f) contributions to E^^'^'^. The optimal error 
estimate taking into account all of these contributions is obtained by choosing 
the scaling r ~ ^/r, which leads to an overall error of order 0{^/r). 
As already observed, formula (5) predicts bending modulus 

= (6) 

for any geometry, and Gaussian rigidity = —^kh/'^, differently from 
what found in Seung and Nelson, 1988 and Gompper and Kroll, 1996, and 
in agreement with the estimates of kh and given in Lidmar et al., 2003. 

5 Continuum model and elastic moduli 

Let us introduce the total limiting energy density 

iy"'"''(I(u),II(n)) = W"'''^{I{u)) + iy^'^"'^(I(n),II(n)) 
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where in the regime of finite elasticity with hnear springs, i.e., W{s) = — 

yj^stretch , 



I(n)) 



^S^fii + 35-22 + 6 - 2\J gii- 2V^gi2 + 85-22 

- 911 + ^VSgu + 3^22 - 4:y/g^^ , 

= ^o( 2, iaiiihl, + 2hl, - 2huh22 + S/i^J 

12(511^22 - 512) 

-8512/111/112 + 2g22{h\^ + 3/112))- 



By differentiating I^'™^* with respect to g = (511, 522, 2512) and h, we get 
the following expressions for the referential stretching stresses N; bending 
stresses M; stretching moduli A; and bending-stretching moduli B 



N 



dW 



limit 



9g 



M 



dW 



limit 



dh. 



i2 Ti/limit 



4- 



i2 Ti/limit 



B 



2 Ti/limit 



D 



which define the limiting elastic response of the network, in association with 
the referential bending moduli (4). The overall (stretching-bending) elasticity 
matrix consists of the second order 6x6 tensor 



C 



A 


B 


B^ 


D 



(7) 



We provide in an Appendix the expressions of the arrays N, M, A, B for 
the general case of finite elasticity, which corresponds to fully anisotropic 
response of the network (full rank C, cf. also Dischcr ct al., 1997). 

Let us introduce now the inverse C"^ = {C~j^) of C, and the 'engineering' 
elastic coefficients 



1 



+ 2Ci2 + C; 



-1 ' 



G 



1 



12 



22 



-0-1 ' 



(8) 



C44 + C55 + 2C45 + C*66 

= : 



V3D (511 + 2522) 
6detl 



(9) 
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representing the area compression modulus, in-plane shear modulus and bend- 
ing modulus of the hmiting surface, respectively. The definitions (8) are con- 
sistent with the notation usually employed for the in-plane engineering con- 
stants of anisotropic plates (refer to standard textbooks on the mechanics of 
composite materials, such as, e.g.. Nettles, 1994; Sec. V.C). For what instead 
concerns the bending modulus (or rigidity) kh, we recall that the theory of 
curvature elasticity of fluid membranes defines such a quantity as the rigidity 
associated with the sum of the principal curvatures, assuming isotropic bend- 
ing response (Hclfrich, 1973). Within the present anisotropic framework, it 
is easily verified that (9) conventionally defines kh as the arithmetic mean of 
the bending rigidities to the incremental (virtual) strains with g' = g" = 0; 
h' = {1/2,1/2,1} and h" = {1/2,1/2,-1} from the current configuration 
(noticeable strain rates that modify the trace H = Ti II = hi -\- h2 of the cur- 
vature tensor, leaving the Gaussian curvature K = detll = hih2 — {h^/2)'^ 
and the first fundamental forms unchanged). We wish again to remark that 
the elasticity matrix (7) depends both on the first and second fundamental 
forms of the current (deformed) configuration of the limiting surface (cf. the 
formulae given in the Appendix), introducing a type of stretching-bending 
coupling that is not accounted for in standard energy functionals of fiuid 
membranes (refer e.g. to Lipowski and Girardet, 1990; Helfrich and Kozlov, 
1993). 

We now consider the bending- dominant regime with I ^ Id. It is not 
difficult to verify that, in such a regime, equations (8) reduce to 

I^A = -TT- (10) 



lQV3D{hli + 3/i?2) + 12\/3fc ^ 

- 2hl^h22 - 6hiihl^h22 + 9hl^hl^ - 3hl^{hl^ - hl2)) 

+ 6Dk {3hl-^ + Shl^ - 2/iii/i22 + ^hl^) + 9^^) 



Formula (10) is in agreement with the result presented in Zhou and Joos 
(1997), Scung and Nelson (1988), Disclicr ct al. (1997) for the area compres- 
sion modulus of triangular nets in the small strain regime. On the contrary, 
formula (11) differs from the analogous one provided in the same works for 
the in-plane shear modulus in the regime under examination (/i = y/Sk/A, 
assuming pure planar deformation), due to the presence of curvature terms. 
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The presence of such terms in the expression of Gu imphes that the shear 
response of the examined membrane network is actually anisotropic and cur- 
vature dependent, when marked bending deformation occurs. In the fully 
infinitesimal deformation regime (I ~ Id and Z}|hp <^ k) we instead recover 
the result /i = G12 = y/Sk/A from (11). It is also immediate to verify that 
formula (9) reduces to (6) when gu ^ (722 ~ 1. 



6 Special cases 

We examine in the present section the two special cases of cylindrical and 
spherical networks covered with A^a equilateral triangles. In both such cases, 
we show that formula (5) exactly reproduces the limit value of the total 
dihedral energy of the network for Na 00. The given examples, which 
allow for a comparison with explicit formulae in the literature, are aimed 
to analytically illustrate the universal nature of formula (5), which predicts 
equal bending stiffnesses for any geometry, namely kh = \/3D/2 and = 
—Akh/'^ = —2\/3D/3, at variance with the prediction kh = —hg = 
of Seung and Nelson (1988). We also provide the expressions of the in-plane 
shear moduli of the same networks accounting for finite curvature. The 
area compression modulus is given by equation (10) in the regime under 
consideration, for any geometry. 



6.1 Cylindrical network 

Gomppcr and Kroll (1996) compute the following value of the dihedral en- 
ergy per unit of length of an infinitely long cylinder of unit radius covered 
with equilateral triangles (again per unit of length) 



^f;E-r^4^.(i-j^+o(i/ivi)) 



which reduces to (5) with = 1, K = and \U\ = 27r. Examining now the 
in-plane shear modulus, and considering a cylinder of finite radius R, from 
(11) we deduce the result 

_ VSk V3DH^ _ V3k VSD 
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6.2 Spherical network 

The total dihedral energy of a spherical network covered with TVa equilateral 
triangles is given by formula (19) of Gompper and Kroll (1996): 

A-kD ( , 27r „ ,,.ox\ 47rL' 



73' 



This exactly coincides with (5) for the sphere with H = 2 / K = 1 / and 
\U\ = AnR"^. For what concerns the in-plane shear modulus of a spherical 
network, from (11) we easily deduce 

A- - G,. = — + = _ + _. (12) 

7 Experimental characterization of material 
constants 

Thinking of cell membranes, and especially of blood cell membranes, the re- 
sults presented in the previous sections are useful to conduct an experimental 
characterization of the material constants k and D through methods like mi- 
cropipette aspiration, optical tweezes and fast phase constrast microscopy, 
among others (refer, e.g., to Mofrad and Kamm (2006), Chapt. 2). 

In particular, micropipette aspiration tests are often used to determine 
the elastic properties of nonadherent cells, through measurements of the cor- 
tical tension and increase in the surface area (area compression modulus, 
cf. Evans et al., 1976; Evans and Waugli, 1977); the length of the cell tongue 
aspirated into the pipette at increasing values of applied pressure (shear mod- 
ulus, cf. Waugh and Evans, 1979); the critical aspiration pressure at which 
buckling of aspirated cell tongue occurs (bending modulus, cf. Evans, 1983); 
and other different quantities (refer to the above mentioned works and therein 
references). Two measurements are sufficient to determine k and D (typically 
the area compression modulus and/or the shear modulus, plus the bending 
modulus), while a third one can be useful for a cross-check. Fig. 2 illustrates 
the aspiration of a red blood cell into a micropipette. When the size L of the 
aspirated cell tongue is such that this portion of the cell assumes a spherical 
shape (Fig. 2), and contemporarly stretching strains are sufficiently small. 
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one can predict the shear modulus through (12), with R ^ L. When, instead, 
the curvature of such an element is associated with large stretching strains, 
one should numerically predict the elastic moduli through equations (8) and 
(9), making use of suitable measurements of the stretching strains. 




Figure 2: Aspiration of a red blood cell into a micropipette. 



8 Conclusions 

We have developed universal formulae for the limiting elastic energy of static 
membrane networks, which arise in the limit for the number of vertices/particles 
tending to infinity. The given formulae account for arbitrary geometries, 
large deformations and bending-stretching coupling. Their particularization 
to the infinitesimal elastic regime has been carried out, and comparisons with 
results available in the literature have been established (Scung and Nelson, 
1988; Gompper and Kroll, 1996; Dischcr ct al, 1997; Zhou and Joos, 1997; 
Lidmar et al., 2003). It has been shown that the formulae here presented 
overcome some deficiencies of earlier results for triangulated networks en- 
dowed with cosine-type dihedral angle potentials, which are related to shape- 
dependence of the bending stiffness and sensitivity to mesh distortion of 
the limiting bending energy (Miiller et al., 2006; Gompper and Kroll, 1996; 
Nelson et al., 2004; Espriu, 1987; Baillie at al., 1990). We remark that the 
finite-elasticity version (3) of the present limiting bending energy is able to 
handle heavily distorted triangulations, due to the presence of finite mem- 
brane deformation terms. We have derived the entire set of the elastic co- 
efficients of the examined network model, obtaining general expressions of 
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such quantities, which depend on the first and second fundamental forms gij 
and hij of the hmiting surface (full bending-stretching coupling). We have 
also obtained simplified expressions of the area compression modulus, in- 
plane shear modulus and bending modulus in the bending-dominant regime, 
which are useful to carry out the experimental identification of such quan- 
tities. We address to future work the derivation of deformation-dependent 
elastic moduli of fluctuating membrane networks undergoing large stretching 
and/or bending deformations, on combining the formulae for the limiting 
energies and elastic moduli here proposed with MD or Monte Carlo simu- 
lations, and the local maximum-entropy approximation scheme presented in 
(Fraternali et al., 2010); an extension of the study of phase transition phe- 
nomena for Hookean spring networks presented in Discher et al. (1997), in 
order to account for the coupling of in-plane and out-of-plane deformations 
in the finite elasticity regime; as well as the experimental characterization 
of the examined model through established laboratory tests, and innovative 
techniques to be designed via computer simulation. 
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Appendix. Supplementary data 

The streching and bending stresses and moduli from Sect. 5 can be easily 
calculated with a computer algebra program. Below we give a mathematica 
code to compute these quantities. A complete list of these constants can 
be found as supplementary data associated with this article in the online 
version. 

Introduce the energy density as: 

WlimG[{gl_, g2_, g3_, hl_, h2_, h3_}] 
= (k / (2 Sqrt[3])) 

* (3gl+3g2+6-2 Sqrt[gl - 2 Sqrt [3] g3 /2 + 3 g2] - 
2 SqrtCgl + 2 Sqrt [3] g3/2 + 3 g2] - 4 Sqrt[gl]) 
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+ ((Sqrt[3] D) / (12 (gl g2 - g3'2/4)) 

* ( gl (hl"2 + 2 h3'2/4 - 2 hi h2 + 3 h2'2) - 8 g3 hi h3/4 + 

2 g2 (hl'2 + 3 h3"2/4) )) ; 

WlimE[{el_, e2_, e3_, hl_, h2_, h3_}] 

:= Wlim[{2 + 2 el, 2 + 2 e2, 2 e3, hi, h2, h3}] ; 

The 6 dimensional vector of the stretching and bending stresses (N, M) 
is then given by: 

D[WlimE[{el, e2, e3, hi, h2, h3}] , {{el, e2, e3, hi, h2, h3}, 1}] 

/. {el ^ 1 - gll/2, e2 ^ 1 - g22/2, e3 -gl2, 

hi ^ hll, h2 ^ h22, h3 ^ 2 hl2} // MatrixForm // Simplify 

The 6x6 stretching- bending elasticity matrix C from (7) is given by: 

D[WlimE[{el, e2, e3, hi, h2, h3}] , {{el, e2, e3, hi, h2, h3} , 2}] 

/. {el ^ 1 - gll/2, e2 ^ 1 - g22/2, e3 -gl2, 

hi ^ hll, h2 ^ h22, h3 ^ 2 hl2} // MatrixForm // Simplify 
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